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ABSTRACT 

In this work the thermal emission over cm to sub-mm wavelengths from the winds 
in short-period O+O-star binaries is investigated (potential non-thermal emission is 
presently ignored). The calculations are based on three-dimensional hydrodynamical 
models which incorporate gravity, the driving of the winds, orbital motion of the stars, 
and radiative cooling of the shocked plasma. The thermal emission arises from the 
stellar winds and the region where they collide. We investigate the flux and spectrum 
from a variety of models as a function of orbital phase and orientation of the observer, 
and compare to the single star case. The emission from the wind-wind collision region 
(WCR) is strongly dependent on its density and temperature, being optically thick 
in radiative systems, and optically thin in adiabatic systems. The flux from systems 
where the WCR is highly radiative, as investigated for the first time in this work, can 
be over an order of magnitude greater than the combined flux from identically typed 
single stars. This excess occurs over a broad range of wavelengths from cm to sub-mm. 
In contrast, when the WCR is largely adiabatic, a significant excess in the thermal 
flux occurs only below 100 GHz. 

In circular systems with (near) identical stars the observed variability in synthetic 
lightcurves is typically less than a factor of 2. Eccentric systems may show order of 
magnitude or greater flux variability, especially if the plasma in the WCR transitions 
from an adiabatic to a radiative state and vice- versa - in such cases the flux can display 
significant hysteresis with stellar separation. We further demonstrate that clumping 
can affect the variability of radio lightcurves. 

We investigate the spectral index of the emission, and often find indices steeper 
than +0.6. Synthetic images display a variety of morphologies, with the emission 
sometimes resembling an intertwined "double-helix" . We conclude by comparing our 
results to observations. The predictions made in this paper are of interest to future 
observations with the next generation of radio and sub-mm telescopes, including the 
EVLA, e-MERLIN, ALMA, and the SKA, and future upgrades to the VLBA. 

Key words: shock waves - stars: binaries: general - stars: early-type - stars: mass 
loss - stars: winds, outflows - radio continuum: stars 



1 INTRODUCTION 

Luminous stars of OB and Wolf-Rayet (WR) spectral types 
drive fast, dense, ionized winds with mass-loss rates that are 
high enough to significantly affect their evolution and their 
surroundings. The winds emit free-free radiation over a wide 
range of wavelengths, from cm emission in the radio band, to 
/im emission in the near-infra-red. Since the bremsstrahlung 
opacity is proportional to the square of the wavelength, 
the observed emission originates above a characteristic 
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radius that increases with wavelength (typical values are 
~ 50 - 100 R, at 20 cm (1.5 GHz), and ~ 2 - 5R* at 1mm 
(300 GHz)). At wavelengths where the wind is optically 
thick, the thermal emission can be used to determine stellar 



(e.g. lAbbott et al. 
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200 ll; ISchnerr et all 


2007f), since the wind can be assumed 



to be at terminal velocity, and the interpretation of the 
flux is not strongly dependent on details of the ionization 
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conditions (|Wright fc Barlowlll975l ; IPanagia fc Fellil Il975h . 
Wind asphericity typically does not alter radio derived 
mass -loss rates by more than a factor of two l|Schmid- Burgkl 
1 19821 ). 

The free-free flux is expected to behave like 
S v oc u a , where th e spectral index a ~ +0.6 
at cm wavelengths (|Wright fc Barlowl Il975l ) . How- 
ever, the spectral index of "thermal" source^] at 
shorter (mm and far-infrared) wavelengths in c reases 
to around +0.7 t o +0 .9 ('e.g. IWilliams et al.l Il990l: 
Leitherer fc Robertl Il99ll; lAltenhoff, Thum fc Wendkerl 



19941 : iNugis. Crowther fc Willis! Il998l ) Several possible 
causes for this steepening have been suggested, including 
a decrease of the wind temperature or the ionization 
state wi th radius or a n incre a se of the wind spe e d wit h 
radius l)Schmid- Burgkl 1 19821 : iLeitherer fc Robertl Il99ll ). 
However, since temperature gradients only weakly affect 
the flux through the gaunt factor, while the wind speed 
changes relatively little between the cm and mm emitting 
regions, these effects are unable to explain the difference 
between theoretical and o bserved spectral indi ces between 
cm and mm wavelengths. ISchmid-Burgkl |jj)S2 ' also noted 
that deviations from the expected behaviour "should not 
easily be blam e d" on non-s pherical winds. Thi s led b oth 
IWilliams et~all (|l990l ) and ILeitherer fc Robertl (|l99ll ) to 
suggest a gradient in ionization as the most likely cause. 

Another possibility involves a radial variation of the de- 
gree of clumping within the wind. For instance, an excess 
of mm emission over that expected from scaling the radio 
emission to shorter wavelengths can be caused by a decline in 
the degree of clumping between t he regions of the wind emit- 
ting at mm and cm wavelengths dRunacres fc Blomme|[T996l : 
iBlomme et"aH 120021 ). INugis et al.l (| 19981 ) have also argued 
for models which combine clumping and a highly ionized 
outer wind region. A recent review o f small-scale structure 
in massive star winds can be found in IPuls. Vink fc Naiarrol 
(|2008l ). 

At very high frequencies, the photospheric emission, 
which is rapidly rising, begins to dominate the emission from 
the wind. This transition occurs at roughly the same fre- 
quency as the whole wind becomes optically thin. The in- 
trinsic thermal spectrum from the wind only then flattens, 
as the emission then arises from essentially a blackbody of 
constant projected linear radius R c at all f requencies greater 
than the critical value given by Eq. 12 of I Wright fc Barlowl 
(|l975l ). However, the dominance of the photospheric emis- 
sion means that observationally the emission retains a strong 
wavelength dependence. 

The presence of a companion star may also change the 
observed thermal flux, its spectral index, and the inferred 
mass-loss rates. First, there are now two stars, rather than 



1 We repeat that in this work we concentrate solely on the ther- 
mal emission from the winds of early-type stars, though obser- 
vations reveal that there is also a significant population of non- 
thermal emitters which appears to be linked to the pr esence of 



a companion star l lDouehertv &: Williams 2000: |Pe Beckei 
and the resulting collision of the winds (see |PittarcT^r*al 



2007) 



2001 



and references therein, for models of the non-thermal emission). 
The non-t hermal emission from the hydrodynamical models pre- 
sented in iPittardl d2009al) , which are the basis for the current 
work, will be investigated in future papers. 



one, and the winds of both will contribute to the observed 
emission. Moreover, the interaction o f the winds m ay also 
produce significant free-free emission. IStevend (|l995l ) inves- 
tigated the total thermal flux from a wide colliding winds 
binary (CWB) as a function of the wind momentum ratio, 
and discovered that it could be about 50 per cent higher than 
that expecte d from the m o re ma ssive wind alone. In more 
recent work, IPittard et al.l (|2006l ) investigated the emission 
from an adiabatic WCR, where the plasma in the WCR re- 
mains at high temperature as it flows out of the system, as 
a function of the stellar separation, d acp . If the WCR re- 
mains optically thin, the thermal flux scales as d~p. More- 
over, since optically thin emission scales as S v oc u~° , the 
flux from the WCR may dominate the thermal flux from the 
system, particularly at cm wavelength^. However, there has 
been no investigation to date of the thermal emission from 
systems where the WCR is highly radiative. This situation 
clearly needs to be addressed. 

In this paper we investigate the thermal radio to sub- 
mm emission from short period O+O-star binaries com- 
puted from the three-dimen sional hydrodynamical models 
described in IPittardl (|2009al . hereafter Paper I). We aim to 
illuminate the relative contribution of the WCR to the total 
free-free flux, the effect that this may have on the observed 
spectral index, and the size and nature of orbital related 
variations, for a variety of models where the behaviour of 
the plasma in the WCR is either highly radiative or largely 
adiabatic. 

The layout of this paper is as follows. We briefly de- 
scribe in Section [2] the hydrodynamical models and in Sec- 
tion 13 our method of calculating the thermal cm to sub-mm 
emission. In Section [4] we present our results, focussing first 
on calculations of the emission from the wind of a single 
star, and then on calculations based on our binary models. 
In Section [5] we compare our findings against observations, 
and in Section [6] we summarize and conclude this work. 



2 SUMMARY OF THE HYDRODYNAMICAL 
MODELS 

The radio calculations in this paper are based on the three- 
dimensional hydrodynamical models described in Paper I. 
The models incorporate the radiative driving of the stellar 
winds, gravity, orbital effects, and cooling, and are summa- 
rized in Tables [1] and [2] The models were not designed to 
simulate particular systems, but were aimed at mimicking 
much of the interesting phenomena that can occur. 

Model cwbl is of an 06V+06V system with a circular 
orbit and a period of 3 days. The WCR is highly radia- 
tive, and significantly distorted by orbital effects, showing 
strong aberration and do wnstream curvatur e. Model cwbl 

12007. and references 



is similar to DHCep 
therein), HD 165052 



Linder et al 



Arias et al 



and HD 159176 l|De Becker et al 



2002 



2004bl : 



Linder et ail 1 20071) . 



Linder et aLll2007l ). 



all of which have near identical main-sequence stars of spec- 
tral type O6— 07, and circular or near-circular orbits with 



2 Since this emission has a negative spectral index i t can mimic a 
synch rotron emission component in the system - see lPittard et ahl 



||2006|) for further details. 
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Table 1. Assumed binary parameters for the models calculated in Paper I. The semi-major axis is 34.26 Rq in model cwbl, 76.3 R,0 
in models cwb2 and cwb3, and 55 R© in model cwb4. e is the orbital eccentricity and r) is the (terminal velocity) momentum ratio of 
the winds. v OI ^ and u w are the orbital speeds of the stars and the preshock wind speeds along the line of centres. \ is the ratio of the 
cooling time to the characteristic flow time of the hot shocked plasma. \ 1 indicates that the shocked gas rapidly cools, while X ^ 1 
indicates that the plasma in the WCR remains hot as it flows out of the system. Larger values of the ratio !) OI |,/ii w produce a greater 
aberration angle, ab , and tighter downstream curvature, of the WCR. The degree of downstream curvature of the WCR in the orbital 
plane is given by « cor i ii a , where the curvature is assumed to trace an Archimedean spiral which in polar coordinates is described by 
r = Q coriolis$- The value of a cor iolis corresponds to the approximate downstream distance (in units of d Bep ) along the WCR for each 
radian of arc it sweeps out in the orbital plane. Smaller values indicate tighter curvature. The leading and trailing arms of the WCR 
in model cwb3 display differing degrees of curvature, so the value quoted for this model is an average. The pre-shock orbital and wind 
speeds in model cwb3 are also different for each star/wind - the first (second) value is for the primary (secondary) star/wind. The values 
of X, v OI b/Vv,, a b and Ocoriolis are phase dependent in model cwb4, because of its eccentric orbit - values at periastron and apastron 
are quoted. The values for a cor iolis are calculated after comparing the orbital speeds at periastron and apastron against those in models 
cwbl and cwb2, and represent the "instantaneous" curvature at these phases. 
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0.40 


17 
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cwb2 


O6V+O6V 


10 


0.0 


1 
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1630 
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0.14 


3-4 
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cwb3 


O6V+O8V 
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0.0 


0.4 
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0.16 - 0.084 


~ 2 


4.5 


cwb4 


O6V+O6V 


6.1 


0.36 


1 


334 - 156 


710 - 1665 


0.34 - 19 


0.47 - 0.09 


21-4 


3-10 



periods near 3 days. Hence the hydrodynamics of, and emis- 
sion from, the WCR in model cwbl should be a reasonable 
approximation to these systems. 

Model cwb2 has identical parameters to model cwbl 
except for an increase in the orbital period to 10 days. The 
winds now have more room to accelerate before they col- 
lide, and the postshock gas remains largely adiabatic as it 
flows out of the system. Both the aberration angle and the 
downstream curvature of the WCR are lessened relative to 
model cwbl. Model cwb2 is similar to HD93161A, an 08V 
+ 09V syste m with a circular orbit and an orbital period 
of 8.566 days jNaze et al.ll2005l ). albeit with slightly more 
massive stars and powerful winds. Another sys tem which is 
not too diss imilar is Plaskett's star (HP 47129. iLinder et alj 
I2006J . I 2008), though this object contains stars which have 
evolved off the main sequence. 

Model cwb3 examines the interaction of unequal winds 
in a hypothetical O6V + O8V binary. The stars have the 
same separation as those in model cwb2, but a slightly longer 
orbital period. The stronger wind from the primary star col- 
lides at higher speeds than the slower secondary wind, result- 
ing in postshock plasma which is hotter and more adiabatic. 

Model cwb4 investigates the effect of an eccentric orbit, 
which takes the stars through a separation of 34.26—76.3 R0 
(i.e. the separations of the stars in the circular orbits of mod- 
els cwbl and cwb2). The WCR is radiative at periastron 
and adiabatic at apastron, and its aberration and down- 
stream curvature are phase dependent. A surprising find- 
ing from Paper I is that dense cold clumps formed in the 
WCR at periastron still persist near the apex of the WCR 
at apastron. This is because the clumps have relatively high 
inertia, and require a significant fraction of the orbtial pe- 
riod to be accelerated out of the system. Some well-known 
O+O binaries with eccentric orbits include (in order of in- 
creas ing orbital period) HP 152248 (e = 0-127: ISana et all 
|2004| ), HP 9 3205 (e = 0.46: iMorreU et alj 120011 ). HP 93403 
(e = 0.234; iRauw et all 120021), CvgOB2#8A (e = 0-24 ; 

Pe Becker et all 120061 ) 

e ■ 



Table 2. Assumed stellar parameters for the models. 



|Pe Becker. Rauw fc Manfroidl 12004: 
and tOrionis (e = 0.764: iBagnuolo et alj|200ll ) 



Parameter /Star 


O6V 


O8V 


Mass ( Mq) 


30 


22 


Radius ( R0) 


10 


8.5 


Effective temperature (K) 


38000 


34000 


Mass-loss rate ( Mq yr" -1 ) 


2 x KT 7 


10- 7 


Terminal wind speed (kms - 1 ) 


2500 


2000 



MODELLING THE RADIO EMISSION 
FROM CWBS 



We use the model d eveloped bv lPoughertv et~aL I (|2003l) and 
iPittard et al.l (|2006l ) to determine the thermal radio to sub- 
mm emission from each of our colliding winds simulations. 
Petailed notes on our method can be found in these papers, 
but the salient points are summarized below. 

To calculate the radio emission we read our hydrody- 
namical models into a radiative transfer ray-tracing code, 
and calculate appropriate emission and absorption coeffi- 
cients for each cell. A synthetic image on the plane of the sky 
is then generated by solving the radiative transfer equation 
along suitable lines of sight through the grid. The thermal 
emission (eJJ) and absorption (a^) coeffici ents at frequency 
v are given bv lRvbicki fc Lightmanl (1 19791 ) as 



6.8 x 10~ 38 Z 2 n e rnT~ 



1/2 -hv/kT 

e g s , 



*5 = 3.7 x 10 8 T~ 1/2 Z 2 ?i c ni^ 3 (l - e 



hv/kT\ 



(1) 

(2) 



where Z is the mean ionic charge, n e and m are electron and 
ion number densities, T is the tempera ture of the gas , and 
gs is a velocity averaged Gaunt factor (|Hummerlll988l ). All 
quantities are in cgs units. The densities n e and ni are de- 
termined from the cell density, composition and ionization, 
where the ionization is specified as a function of cell temper- 
ature. The appropriate ionization for a particular tempera- 
ture regime is determined and the absorption and emission 
coefficients are then evaluated. The unshocked winds are as- 
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Table 3. The characteristic radius of emission, R c , and thermal 
flux, S„, of the wind of the 06V star as a function of frequency 
(see IWrieht fc Barlowl Il975l) , assuming that H and He are both 
singly ionized, that the wind temperature is 10000 K, that the 
wind is smooth (i.e. no clumping), and that the star is at a dis- 
tance of 1 kpc. The relatively low values of R c and S v reflect 
the adopted mass-loss rates of the stars (which arc nevertheless 
consistent with recent downward revisions of the actual mass-loss 
rates of massive stars), coupled with our assumption of smooth 
winds. Assuming v(r) = Voo reduces R c and Sv below their actual 
values, particularly at the highest frequencies which probe the in- 
nermost parts of the acceleration region (in fact , the formal values 
of R c at 1 and 2THz obtained from Eq. 11 of I Wright fc Barlowl 
(|l975h are < R»). For the 08V star, R c and S i/ are smaller by 
factors of 0.731 and 0.534, respectively. 



Frequency 


Wavelength 


Rc 


s v 


(GHz) 


H 


(R©) 


(mjy) 


5 


60000 


386 


0.0049 


8.4 


35714 


273 


0.0067 


15 


20000 


185 


0.0094 


22.5 


13333 


111 


0.012 


43.3 


6928 


92 


0.017 


100 


3000 


44.1 


0.028 


150 


2000 


32.9 


0.035 


250 


1200 


22.7 


0.046 


500 


600 


13.6 


0.067 


1000 


300 


8.13 


0.095 


2000 


150 


4.80 


0.13 



sumed to be isothermal at T = 10000 K, and each model is 
assumed to be at a distance of 1 kpc. 

In the majority of the models in this work we assume 
that the winds are smooth (i.e. not clumped). This enables 
the characteristic radius of emission and the r = 1 surface to 
better fit within our hydrodynamical grids, at the expense 
of underestimating the emission compared to real systems 
where the winds are often strongly clumped. Our assump- 
tion of smooth winds also removes the need to specify the 
radial stratification of the clump vo lume filling facto r, /, 
which remains poorly known (see e.g. IPuls et al.|[2008l ). Ta- 
ble |3] lists the characteristic radius, R c , and flux density, 
S v , of free- free emission from our model 06V star at vari- 
ous frequencies between 5 and 2000 GHz (corresponding to 
wavelengths between 6 cm and 150 fj.m) . The radius of opti- 
cal depth unity, f? T =i = 0.623 R c - 

All of the hydrodynamical grids in our models are large 
enough to contain the characteristic radius for free-free emis- 
sion at 43 GHz, while models cwb2 and cwb3 are also able 
to contain the characteristic radius for emission at 15 GHz 
(compare Table [3] with Table 3 in Paper I). This enables us 
to explore the thermal radio emission from these systems, in 
addition to the X-ray emission explored in a previous paper 
(|Pittardll2009bl . hereafter Paper III). 

Clumping increases the emission and absorption coef- 
ficients in Eqs. [T] and [2] by a factor of 1//, where / is the 
volume filling factor of the clumps, and increases the char- 
acteristic radius by a factor / _1 ^ 3 , and the flux density by 
a factor / _2//3 . This means that the flux from a model with 
a given M and /, if interpreted as arising from a wind with 
no clumping, gives (incorrectly) the mass- loss rate M /yff 
(note that one can also take an alternative approach: a par- 



ticular observed flux can be translated into a mass-loss rate 
Afsmooth assuming that there is no clumping, or a mass-loss 
rate M c i umpo( j = M smoo thV7i assuming that there is clump- 
ing). However, even when the stellar winds are clumpy, the 
WCR may be considerably less so. This is expected to be the 
case when the postshock gas behaves largely adiabatically. 
In such situations the destruction time of the clumps can 
be significantly shorter than the flow time of the shocked 
gas out of the system, and the WCR, though highly tur- 
bulent, is no ticeably smoother than the s tellar winds (see 
Pittard||2007l). Therefore, i n previous work (iDoughertv et al.l 



20031 ; I Pittard et al. | 120061 ; iPittard fc Dougherty! 120061 ) we 



have assumed that the winds are clumpy and the WCR is 
smooth(er). 

In this work we also investigate the emission assuming 
clumpy winds with / = 0.1, but a WCR which is smoother 
to varying degrees i.e. 0.1 ^ /wcr ^ 1. This provides more 
realistic fluxes from the unshocked winds and insight into 
how clumping in different regions of the system affects the 
observed lightcurves and spectra. In these models the coef- 
ficients in Eqs. [1] and [2] are multiplied by a factor of l// w inds 
(1//wcr) for cells in the unshocked winds (WCR), respec- 
tively. 



4 RESULTS 

4.1 Emission from the wind of a single star 

We first investigate the cm to infrared free-free emission 
from the wind of a single star, namely from the 06V star 
used in our binary simulations. The wind density as a func- 
tion of radial distance is o btained from our solution to 
the modified CAK equ ations (|Castor. Abbott fc Kleirill 19751 : 
iPauld rach et all Il986l ). At high frequencies the emission 
arises from within the acceleration zone o f the win d, and it 
is no longer possible to use, e.g., Eq. 8 in I Wright" 



Barlowl 

(|l975l ) to evaluate the flux. We therefore numerically inte- 
grate the flux as a function of impact p arameter, follow- 
ing the method of IWright fc Barlowf l| 19751 ). The wind is as- 
sumed to begin at the point where the outflow velocity is 
equal to the isothermal sound speed, which is w 25 kms -1 
in our models. As in the binary models which follow we as- 
sume that the wind is smooth, isothermal and maintains a 
constant ionization: these assumptions allow us to isolate 
the effect of the wind acceleration on the resulting emission. 
We also include the photospheric flux, which is modelled as 
a black body for simplicity. The photospheric flux is atten- 
uated at wavelengths where the wind is optically thick. 

Fig. OJa) shows the resulting spectrum obtained from 
the wind and photosphere of our model 06V star. At GHz 
frequencies (cm wavelengths) the total emission is domi- 
nated by the wind and arises from large distances where 
the wind has very nearly reached its terminal speed. The 
spectral index (Fig. Hp) is near +0.6, as expected. The emis- 
sion starts to probe the acceleration region of the wind at 
v > 100 GHz, where a steepening of the spectrum becomes 
apparent. It is from this point that the flux increasingly ex- 
ceeds that from a calculation where the wind is forced to 
emanate from the star at its terminal speed (i.e. instanta- 
neous acceleration). As shown in Fig. [Ha), such a termi- 
nal speed wind becomes optically thin at frequencies above 
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Figure 1. (a) The spectrum from the wind and photosphere of a simple model of a single O6V star. The photosphere is assumed to 
emit as a blackbody at a temperature of 38000 K, while the wind is assumed to be at a constant temperature of 10000 K and to have 
a constant ionization structure of H + , 50 per cent neutral He and 50 per cent He+, and CNO + . The photospheric emission is strongly 
absorbed by the wind for v < 2 X 10 12 Hz. (b) The spectral index of the total and wind component of the spectrum, (c) The variation 
of the wind acceleration index, e, where v(r) oc r e . 




Figure 2. Intensity images from model cwbl at i = 0°. The maximum intensity (black in the images) is 4.6 X 
10~ 10 ergem" 2 s" 1 Hz" 1 ster^ 1 at 43GHz, 2.7 X 10~ 9 ergem" 2 s" 1 Hz" 1 stcr" 1 at 100GHz, 7.1 X 10~ 9 ergem" 2 s" 1 Hz" 1 stcr" 1 at 
250GHz, and 1.1 X 10~ 7 ergem" 2 s" 1 Hz -1 ster -1 at 1000GHz. The gray scale spans 4 orders of magnitude between maximum and 
minimum. The major ticks on each axis mark out 0.2 mas. Orbital phases are marked on for an observer in the orbital plane, with phases 
0.0 and 0.5 corresponding to conjunction and phases 0.25 and 0.75 to quadrature. The stars and the WCR rotate anticlockwise. We also 
define an azimuthal viewing angle, </>, which increases anticlockwise from the bottom of the image, so that at phase 0.0 <j> = 0°. <j> = 90° 
corresponds to phase 0.75, and cA = 180° and <j> = 270° correspond to phases 0.5 and 0.25, respectively. 



700 GHz (Eq. 12 of lWright fc Barlow|ll975r i. with declining 
flux thereafter (see also lBertout et al.lll985T ). In contrast, the 
sharp increase in the wind density of the CAK model as the 
photosphere is approached results in a steep rise in the ther- 
mal flux. At frequencies around 2 x 10 12 Hz (A = 150 /im), 
the spectral index of the emission from the wind steepens to 
reach a maximum of ~ +1.1. This steepening is due purely 
to the observed emission occuring from regions where the 
wind is rapidly accelerating. 

The wind also starts to become optically thin at v ~ 
2 x 10 12 Hz (A = 150 /im), and there is a steep rise at higher 
frequencies in the spectral index of the total emission as the 
photospheric emission becomes dominant at v > 4 x 10 12 Hz 
(A = 75 fim; for stars with denser winds, the emission from 
the wind may dominate the photospheric emission until /im 
wavelengths). The maximum spectral index for thermal free- 
free emission is +2.0, the value for a blackbody source where 
the flux Si, oc v 2 , as given by the Rayleight- Jeans approxi- 
mation to the Planck function. This is indeed the case in 
a calculation where the wind and photospheric tempera- 
tures are equal. However, when there is a sharp change be- 
tween the wind and photospheric temperatures, such as in 
our simple model, spectral indices greater than +2.0 can 
be obtained. This is apparent in Fig. [TJb), where the spec- 
tral index actually reaches +2.5 at v w 5 x 10 12 Hz, due to 
the steeply rising photospheric flux as the absorption of the 



overlying wind declines^ The photospheric emission from 
our model O6V star peaks at v w 3 x 10 15 Hz (A = 100 nm), 
and subsequently declines. 

For a smooth, spherical, isothermal wind of constant 
ionization, the spectral index, a, is rel ated to the accelera- 
tion index, e (where v(r) oc r E ), by (e.g. iLeitherer fc Robertl 
[1993) 



4e+ 1.8 
2e + 3 ' 



(3) 



Note that v(r) oc r E is a local approximation (assumed 
valid over the formation region corresponding to a partic- 
ular wavelength). The value of e is not valid over the whole 
distance range, in the same way that the value of (3 in the 
"/3-law" (see below) is. 

To obtain a = +1.1 requires that e = 0.83. This is 
easily achieved in accelerating winds. Fig. [He) shows e as a 
function of radius for the CAK wind profile and standard 
"/3-law" parameterizations (where 



3 In reality, there will not be such a sharp division in tempera- 
ture between the wind and photosphere, and one would expect 
that the maximum value of the spectral index determined from 
observations would be closer to +2.0. Spectral indices greater 
than +2.0 have in fact been previ ously reported from Be stars 
(iDoughertv. Taylor fc Water j|l99lh . but we note that these were 
based on 3 upper limits to the flux, and so are meaningless. 
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Figure 3. Radio lightcurves of the free-free emission from model cwbl at 43 GHz (top), 100 GHz, 250 GHz, and 1000 GHz (bottom) for 
inclination angles i = 30° , 60°, and 90°. The total free-free emission is shown in the left panels, and the contributions from the winds and 
WCR in the middle and right panels, respectively. In all cases, the observer is located along a direction vector specified by <f> = 0°. The 
stars are at conjunction at phases 0.0 and 0.5, and quadrature at phases 0.25 and 0.75. Note that the y-axes (flux scales) are different 
in each panel. 



v(r) 



+ (voo-vo)(l-R*/ry 



(4) 



and vo is the initial wind velocity at the stellar radius, as- 
sumed here to be 25 kms -1 ). The CAK velocity profile 
is reasonably well approximated by a j3 velocity law with 
f3 = 0.8. Winds with slower acceleration have larger values 
of f3, and vice-versa. Winds with large j3 have smaller val- 
ues of e close to the star, but larger values further from the 
star where their winds continue to accelerate. Values of e 
greater than 100 are achieved in the innermost regions of 
strongly accelerating winds. In the CAK model e — 0.83 at 
r/R* - 1 = 0.97. 

The maximum frequency which can be explored in the 
3D hydrodynamical models is set by the necessity of suffi- 
ciently resolving the acceleration region where the character- 
istic radius of emission occurs. Models cwbl and cwb4 have 
a numerical resolution of 0.5 R©. In models cwb2 and cwb3 
the numerical resolution is 1.25 R©. Therefore, it is only 
possible to examine frequencies up to ~ 1 THz, since the 
optical depth for an observer at infinity to r = 1.1 R* is less 
than unity at frequencies greater than this. In the following 
we restrict our investigation to v ^ 1 THz (A ^ 300 fim) . 



4.2 Emission from binary systems 

4.2.1 Model cwbl 

The hydrodynamical grid used to compute model cwbl con- 
tains within its volume the characteristic radius of emission 
from the 06V wind at v ^ 43 GHz. The optical depth unity 
surface, which is slightly smaller than R c , fits completely in- 
side the grid for v > 15 GHz. Therefore, we limit our study 
of the free-free intensity images and lightcurves from model 
cwbl to v 43 GHz, but show spectra down to 8 GHz for 
ease of comparison with our other models. 

Fig. [2] shows intensity images of the free-free emission 
from model cwbl at 43, 100, 250, and 1000 GHz, for an ob- 
server located directly above the orbital plane (i — 0°). It is 
immediately apparent that the WCR dominates the emis- 
sion, even at a frequency as high as 1 THz. The morphol- 
ogy of the emission from the WCR reflects its underlying 
structure and dumpiness, caused by its high susceptibility 
to dynamical instabilities (see Paper I for a detailed discus- 
sion of the hydrodynamics). Near the apex of the WCR the 
projected emission from different inhomogeneities merges to- 
gether, but further downstream individual clumps and bow- 
shocks can be identified. It is clear that a great deal of po- 
tential emission is not captured due to the limited physical 
extent of the hydrodynamical grid used in the model. For 
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Figure 4. Top: free-free radio spectra from model cwbl at viewing angles: a) i = 90°, <f> = 0°, b) i = 90°, <f> = 90°, c) i = 0°, <j> = 0°. In 
each case the contributions of the free-free emission from the unshocked winds and the WCR to the total free-free emission are shown. 
The stars are at conjunction (phase 0.0) in a), and quadrature (phase 0.75) in b). The stars are also at quadrature in c). There is likely 
a significant loss of flux at the lowest frequencies displayed due to the finite size of the numerical grid. Bottom: the spectral index as a 
function of frequency from the spectra in panels a)-c). 




Figure 5. Free-free radio spectra from models cwbl (a), cwb2 (b), and cwb3 (c), at viewing angles i = 90°, <j> = 0° (phase 0.0, pink), 
and i = 90°, <f> — 90° (phase 0.75, cyan) - conjunction and quadrature, respectively. Because of the symmetry in models cwbl and cwb2, 
the spectra at phase 0.0 and 0.75 are virtually identical to those at phase 0.5 and 0.25, respectively. In model cwb5 the 08V star is in 
front at phase 0.0. 



obvious reasons it is impossible to know just how much flux 
is "lost" in this way as dense clumps in the WCR flow out 
through the grid boundaries, but it is clear that this problem 
is worse at lower frequencies where structure in the WCR 
remains optically thick for larger downstream distance^. 
Nevertheless, these images, and the lightcurves and spectra 
which follow, provide important new insight into the emis- 
sion at cm and sub-mm wavelengths from CWBs. 

The spatial scale of the images is 1.12 by 1.12 mas, 
which is too small to be resolved with the EVLA (angu- 
lar resolution of 4 mas at 50 GHz) or ALMA (highest an- 
gular resolution expected to be about 10 mas). In contrast, 

4 One can estimate a rough upper limit for the amount of flux 
lost by taking the flux at 1000 GHz, and extrapolating to lower 
frequencies with a spectral index of +0.6 - this reveals that the 
actual flux from the WCR at 43 GHz may be up to an order of 
magnitude higher. 



the VLBA has the necessary spatial resolution, but even 
with all current upgrade plans may struggle with the nec- 
essary sensitivity to obtain images of this system within a 
reasonable integration time (e.g. lOhrs) to avoid smearing 
of the essential features due to orbital motion. In 10 hrs, 
with 256 Mbit/s bandwidth at 43 GHz, the VLBA + GBT + 
phased- VLA obtains a sensitivity of 0.029 mjy/beam. Band- 
width improvements to 4Gbit/s will increase the sensitivity 
by a factor of 4. A signal to noise ratio (SNR) of 5 (the min- 
imum needed to provide imaging constraints) in 10 hrs thus 
requires a flux of 0.0725 mjy/beam. The spatial resolution 
at 43 GHz is ~ 0.3 mas, so each arm of the WCR falls within 
a beam (see Fig. [5J), which is roughly 50 per cent of the to- 
tal flux, or 0.02 mjy. Hence a 10 hr exposure is insufficient 
to obtain a SNR of 5 - instead ~ 130 hrs is required. Having 
said this, it has already been noted that there is a significant 
loss of flux at such frequencies from model cwbl due to the 
finite size of the hydrodynamical grid, so the future detec- 
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tion of such a system with the VLBA is not totally beyond 
the realm of possibility (e.g. if the flux from model cwbl 
at 43 GHz was actually 3.6 x higher). Systems with stronger 
winds (higher mass-loss rates) than assumed for model cwbl 
may further improve the possibility of such imaging. 

Fig. |3]shows the free-free emission lightcurves at 43 GHz 
(top panels) to 1000 GHz (bottom panels), as a function of 
the inclination angle, i, of the system. The left panels show 
the combined emission from the unshocked winds and the 
WCR, while the middle and right panels display the un- 
shocked and shocked components, respectively. The stars are 
aligned along the lines-of-sight at phases 0.0 and 0.5, and are 
at quadrature at phases 0.25 and 0.75. The variations in flux 
are entirely due to changes in the occultation and circum- 
stellar absorption, as the intrinsic emission is constant. Note 
that the presence of the companion star breaks the spherical 
symmetry of the unshocked winds. For example, because of 
the twin radiation fields, the winds accelerate more slowly 
towards each other, and more quickly away from each other. 
There are also regions of lower velocity wind in the "shad- 
ows" behind each star (see Paper I for further details). 

The lightcurves of the WCR emission (right panels in 
Fig. [3]l would display dual symmetry about phases corre- 
sponding to both quadrature (0.25, 0.75) and conjunction (0, 
0.5) if there were no orbital effects on the WCR (i.e. no aber- 
ratio n or coriolis induced curvature - cf. IPittard fc Stevens] 
1 19971 ). because of the equal winds and constant stellar sep- 
aration. Orbital effects break the symmetry about quadra- 
ture, so that in the case of identical winds a double peri- 
odicity in the lightcurve is expected. However, dynamical 
instabilities which form in the WCR can break this symme- 
try too, as revealed by close examination of Fig. [3] 

At the lower frequencies the flux from the unshocked 
winds is typically only a small fraction of the expected com- 
bined flux from the winds (see the horizontal red lines across 
the plots in the middle column of Fig. [3}. For instance, at 
43 GHz the thermal flux from a single smooth (unclumped) 
wind of the same parameters as used in model cwbl is 
0.017 mjy, while the combined unshocked winds contribu- 
tion in Fig. [3jb) is ~ 0.01 mjy. This disparity reduces at 
higher frequencies - at 250 GHz, the two unshocked winds 
contribute « 0.055 mjy, whereas the expected single star 
value is 0.046 mjy - and is eliminated by 1000 GHz. 

The lack of observed flux from the unshocked winds as 
shown in Figs. [3jb) and (e) is likely due to the removal of 
a large fraction of the unshocked material from both of the 
winds by the presence of the WCR. In contrast, at the high- 
est frequencies considered (250 and 1000 GHz), the char- 
acteristic radius of emission moves so close to the surface 
of the stars (particularly at 1000 GHz) that the emission 
from this region is not strongly disturbed by the presence 
of the WCR at larger distances. A small additional effect at 
the lowest frequencies arises from the limited extent of the 
computational grid compared to the characteristic radius of 
emission, which causes an underestimate of the flux. This 
can be seen in the spectra in Fig. [4] where a distinct down- 
turn in the unshocked ( "ff winds" ) flux is visible at the lower 
frequencies investigated. However, at higher frequencies be- 
tween 43 — 250 GHz the removal of unshocked wind material 
by the WCR is the main cause of the "missing" flux. 

At 1000 GHz, the combined flux from the unshocked 
winds is ~ 0.25 mjy, which actually slightly exceeds twice 



the theoretical value of 0.095 mjy for our single star. This 
reflects the fact that the characteristic radius of emission is 
now deep within the acceleration zones of the winds, while 
the theoretical flux is estimated assuming that the entire 
wind is moving at its terminal speed. This excess over the 
"theoretical" value was also seen in the single-star calcula- 
tions presented in the previous section. 

The thermal emission from the unshocked winds is gen- 
erally pretty steady with orbital phase. Slight maxima oc- 
cur near to quadrature, with minima near to conjunction, 
as expected. The 1000 GHz lightcurve is not so smooth in 
appearance - its more "noisy" behaviour is a consequence of 
lines of sight through the WCR becoming optically thin in 
some directions, but not in others. The density of the thin 
layer of cooled postshock gas is w 6 x 10~ 14 gem -3 near the 
apex of the WCR. The linear absorption coefficient of this 
gas at 1000 GHz is « 1.5 x 10~ n cm -1 , while the thickness 
of this layer is typically 1 Rq, giving an optical depth of 
approximately unity. Hence we see the emitting regions of 
the WCR and the unshocked winds near the stellar photo- 
spheres through a semi-opaque, semi-porous WCR. 

The free-free emission from the WCR, on the other 
hand, displays much greater phase variability. It has already 
been noted that the flux from the WCR dominates the emis- 
sion from the unshocked winds. Fig. [3]shows that the WCR 
flux is around 5x greater at 250— 1000 GHz, though the true 
value is likely to be somewhat higher. Future calculations 
with a larger grid are necessary to accurately determine this 
ratio. The lightcurves of the emission from the WCR typi- 
cally display two minima per orbital cycle, with a phase-lead 
with respect to conjunction of the stars which appears to 
increase with frequency. Whereas the X-ray lightcurves of 
model cwbl in Paper III produce reasonably sharp maxima, 
the lightcurves in the left column of Fig. [3] display much 
broader maxima. 

An explanation for the behaviour of the lightcurve of 
the emission from the WCR now follows. Figs. [2] and [3] show 
that at radio to sub-mm frequencies, the most intense emis- 
sion from model cwbl occurs from the cold, dense sheet 
of post-shock material of cooled gas within the WCR. The 
intrinsic emissivity of this region is very high, since the free- 
free emission coefficient etf oc p 2 T~ 1 ^ 2 , and the density in 
the WCR is orders of magnitude greater than in the un- 
shocked winds. The WCR is also optically thick, and, be- 
cause of the downstream curvature induced by the coriolis 
effect, it displays a larger projected surface to the observer, 
and hence a larger flux, when the stars are at quadrature 
than at conjunction (see Fig. [2j . The amplitude of the vari- 
ability decreases slightly with frequency as the optical depth 
through the WCR declines. 

Fig. [4] shows radio-to-sub-mm spectra as a function of 
orbital phase for an observer in the orbital plane (panels a 
and b), and directly above/below it (panel c). In all cases 
the optically thick emission from the WCR dominates the 
observed flux (remember that there is likely a significant 
loss of flux at the lowest frequencies displayed). However, 
the emission from the WCR is less dominant over the un- 
shocked winds for an observer above or below the orbital 
plane (compare panel c to panels a or b), since the orien- 
tation of the layer of cold dense gas to the observer is less 
favourable (which reduces the flux observed from it), while 
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that of the winds is more favourable (both winds are seen 
side-to-side on the sky). 

The spectrum of the total thermal flux shows consider- 
able curvature for an observer in the orbital plane (pan- 
els a and b), but significantly less curvature for an ob- 
server with i = 0°. Fig. |3£d)-(f) shows the frequency de- 
pendence of the spectral index from these spectra. Between 
100 and 250 GHz the spectral index of the unshocked winds, 
Qwinds, is +0.95, +0.95, and +1.00, for viewing angles of 
(i = 90°, <j> = 0°), (i = 90°, (t> = 90°), and (i = 0°), re- 
spectively. These are steeper than the canonical +0.6 for 
optically thick, isothermal, terminal speed winds, reflect- 
ing the fact that the winds in our model are accelerating 
and asymmetrical. a w inds increases towards 1000 GHz as the 
emission probes deeper into the acceleration region of the 
winds. At 1000 GHz, a wind s = +1.09, +1.34, and +1.29 for 
the same orientations. The variation of a w inds with viewing 
angle clearly shows that the presence of the companion star 
has a non-negligible impact on this parameter, with large 
variations in a W mds between conjunction and quadrature. 

The reduction in the spectal index at the highest fre- 
quencies of the emission from the WCR, Qwcr, is an indica- 
tion that parts of the WCR are becoming optically thin. Be- 
tween 100 and 250 GHz, a W CR = +1.51, +1.22, and +1.39 
for the three viewing angles noted above. At 1000 GHz, 
Qwcr = +0.66, +0.64, and +0.89, respectively. For the to- 
tal free-free emission between 100 and 250 GHz, a — +1.29, 
+ 1.00, and +1.39, while at 1000 GHz, a = +0.73, +0.77, 
and +1.00. The less dominant contribution of the WCR at 
i — 0° results in the least variation in a with frequency (see 
Fig.Hd)-(f)). 

Fig. [5ja) shows the variation of the total thermal emis- 
sion for an observer in the orbital plane between conjunction 
and quadrature. The emission at quadrature is dominant 
over the emission at conjunction up to v w 500 GHz. Be- 
tween 500 and 1000 GHz the emission at these phases is al- 
most identical. Note that in this frequency range the maxima 
and minima in the lightcurves no longer occur at quadra- 
ture and conjunction. It is not possible to state whether 
this behaviour continues beyond 1000 GHz, for the reasons 
mentioned earlier in Sec. 14.11 

4.2.2 Model cwb2 

Fig- HD shows intensity images of the free- free emission from 
model cwb2 at 15, 43, 100, 250, and 1000 GHz, for an ob- 
server located directly above the orbital plane (i = 0°). The 
curved shape of the WCR is clearly seen. The WCR domi- 
nates the emission at the lower frequencies (v < 100 GHz), 
while the unshocked winds dominate the emission at higher 
frequencies. At 15 GHz the WCR is very bright in compari- 
son to the winds. The brightest part of the WCR is its apex, 
between the stars. Further downstream the gas at the trail- 
ing edge of the WCR is brighter than at the leading edge, 
reflecting the higher densities there (see Paper I). The an- 
gular scale of all of the images in Figs. 16181 is 2.7 by 2.7 mas, 
which is again too small to be resolved with EVLA/ALMA 
but may be resolved with the VLB A. 

Intensity maps from model cwb2 for an observer in the 
orbital plane are shown in Fig. [7] The features in these im- 
ages depend on the relative orientation and the degree of 
emission/absorption of the winds and the WCR, and dis- 



play significant changes as a function of frequency and or- 
bital phase. At conjunction only the foreground wind is seen, 
since this blocks the emission from the wind of the more dis- 
tant star. The surface brightness of the foreground wind is 
less than that of the WCR at 15 GHz, so there is a "hole" in 
the emission from the WCR at conjunction. In contrast, the 
surface brightness of the unshocked winds is greater than 
that of the WCR when v > 100 GHz. The sharp, curved 
edges to the emission from the WCR near the borders of 
the images is due to limb-brightening. If the hydrodynami- 
cal grid were large enough to contain several spirals of the 
WCR around the stars one would expect to see an expand- 
ing series of progressively fainter, concentric, edges. The un- 
shocked winds emission at 1000 GHz arises from a region 
very close to the surface of the stars. 

Some very interesting features arise in the synthetic 
images as the phase of the orbit advances. From the ob- 
server's point-of-view, the foreground star slides to the right, 
which at 15 — 100 GHz exposes more of the brighter WCR 
to the left by phase 0.125. At higher frequencies the winds 
are brighter than the WCR. At phase 0.25 (quadrature) the 
apex of the WCR is seen side on, with the unshocked winds 
to either side. The brightest part of the WCR is between 
the stars as expected, and is visible in these images even at 
1 THz. At low frequencies stronger absorption of the WCR 
occurs on the right-hand-side of the image due to the greater 
depth of intervening wind material. The right-hand column 
of Fig.[7]shows the system viewed at phase 0.375, halfway be- 
tween quadrature and conjunction. The image now displays 
a prominent double-helix-like structure formed through the 
limb-brightening of the WCR. The vertical curvature of the 
WCR induced by the increasing phase-lag of shocked gas 
with height/depth from the orbital plane is strikingly dis- 
played (see also the warped shock surfaces of the WCR 
shown in Paper I). At low frequencies, this double-helix-like 
structure is partially obscured by the foreground wind (now 
to the left of the image). 

Fig. [8] displays images obtained for an observer with an 
inclination angle i — 30°. The brightest parts of the WCR 
are again those which are limb-brightened, with the overall 
morphology of the WCR being shaped like an "S" . In some 
of the images displayed (e.g. the middle-left image with v = 
100 GHz and = 0°) the arms of the WCR do not trace 
back to a central point, due to the additional emission from 
the winds. This figure also demonstrates that the observed 
shape of the WCR, particularly at the lower frequencies, 
owes much to the location of any foreground absorber (i.e. a 
wind), and is not necessarily a good indicator of the actual 
shape and position of the WCR (a point made already in 
iDoughertv et alj|2003h . 

Lightcurves of the thermal emission from model cwb2 
are shown in Fig. [9] Because of the larger size of the com- 
putational volume, lightcurves at 15 and 22 GHz are also 
displayed. The lightcurves of the total observed emission dif- 
fer distinctly in their shape compared to those from model 
cwbl (left panels of Fig. [3} - those from model cwb2 display 
maxima which are less broad. This reflects the differences 
in the optical depth through the WCR in the two models. 
Fig. [10] illustrates that the emission from the WCR of model 
cwb2 becomes optically thin at relatively low frequencies, in 
contrast to that of model cwbl where it remains optically 
thick to v > 1000 GHz. 
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Figure 6. Intensity images from model cwb2 at i = 0°. The stars and the WCR again rotate anticlockwise with time. In this and 
Figs. [7] and [8] the maximum intensity (black in the images) is 1.6 X 10~ 10 ergem -2 s —1 Hz -1 ster -1 at 15, 43 and 100GHz, 3 X 
10 — 10 erg cm -2 s — 1 Hz -1 ster -1 at 250 GHz, and 3.5 X 10 — 9 erg cm -2 s _1 Hz -1 ster" 1 at 1000 GHz. The gray scale spans 4 dex between 
maximum and minimum. The major ticks on each axis mark out 0.5 mas. Some numerical artefacts from the hydrodynamics are visible, 
but have negligible impact on the flux. 



The free-free emission from the unshocked winds at 
43 GHz is still substantially below the combined emission 
from the winds of two single stars (panel h in Fig. [9]l , but 
compared to model cwbl (see Fig. [3]d) the flux from model 
cwb2 is roughly twice as high. This reflects both the larger 
computational domain in model cwb2, which means more of 
the emission at 43 GHz is captured, plus the wider separa- 
tion of the stars i.e. the WCR removes a smaller segment 
of the 43 GHz emitting region from the unshocked winds. 
Once again, it is only by 250 GHz that the flux from the un- 
shocked winds is comparable to twice the flux from a single 
star. And, as was also the case for model cwbl, by 1000 GHz 
the unshocked wind flux exceeds twice the single star value 
(due to the emission probing deep into the acceleration re- 
gion of the winds). However, the 1000 GHz lightcurve of un- 
shocked wind emission from model cwb2 (Fig. (9jj) is almost 
completely constant (the main exception being a triangular 
shaped minimum centered on each conjunction), in stark 
contrast to the greater variability seen from model cwbl 
(Fig. [5]<). The minimum in Fig. (Hq) fails to approach 50 
per cent of the unshocked wind flux. This is expected, since 
two perfectly symmetric winds around stars in conjunction 
will not completely occult, because wind emission from the 
background star will not be completely absorbed for im- 
pact parameters where the winds start to get optically thin. 
Therefore, a perfect occultation by two stars with stellar 
winds is not possible, and hence a dip of 50 per cent is 
never reached. An additional (smaller) part may be played 
by the fact that the winds are no longer spherically symmet- 
ric, since their acceleration is inhibited towards each other, 
but is boosted in other directions. A slight lag in the phase of 
the minima relative to conjunction is increasingly apparent 
at lower frequencies. 

Now consider the observed emission arising from the 
WCR (see the right column of Fig. 0. A maximum in the 
emission again occurs near quadrature, just as in model 
cwbl, though this time for a totally different reason. In 
model cwb2, maxima occur when the line of sight of the 
observer to the apex of the WCR (where the intrinsic emis- 
sivity is greatest) is through one of the low-opacity arms of 
the WCR. This contrasts with model cwbl, where maxima 
occur when the optically thick WCR is oriented almost face 
on. 

Although the free- free emissivity per unit volume in the 
WCR is much smaller in model cwb2 than in cwbl due to 
the reduced densities and higher temperatures, the optical 




250 




Figure 11. Intensity images from model cwb3. The left 
panel shows the emission at 43 GHz with i = 0°. The right 
panel displays the emission at 250 GHz with i = 30° and 
4> = 45° (phase 0.875). The maximum intensity is 1.6 X 
10 -10 erg cm -2 s _1 Hz -1 ster -1 and the major ticks on each axis 
mark out 0.5 mas. = 0° corresponds to the 08V star in front. 



depth within this material is also reduced. The net result 
is that the observable flux from the WCR at v < 50 GHz 
is much higher in model cwb2 than in model cwbl. But 
above 50 GHz, when the opacity through the circumstellar 
environment is low, the reduced densities in the WCR limit 
the intrinsic emission compared to model cwbl, so that the 
observed flux from the WCR in model cwbl exceeds that 
from model cwb2 (compare Figs. l4l and I10|l . 

At 43 GHz the free-free flux from the WCR varies by a 
factor of 2 for observers with a 90° inclination (Fig. [9j). At 
lower frequencies the variation is even more pronounced, as 
the relative contribution of free-free emission from the WCR 
becomes increasingly dominant (the free-free flux from the 
winds scales as u +0 ' (> , whereas it scales as v~°' from the 
optically thin WCR). 

Interestingly, for v ^ 43 GHz, the cwb2 WCR 
lightcurves display a remarkable resemblance to the 0.1 — 
0.5 keV X-ray lightcurves from this model shown in Pa- 
per III, including the asymmetrical slope of the minimum 
and the faster rise to maximum than the fall from it. More- 
over, the 100 and 250 GHz WCR lightcurves (panels 1 and o 
in Fig. |9]) are similar to the 0.5 — 2.5 keV X-ray lightcurve. On 
the other hand, the 2.5—10 keV X-ray lightcurve is most sim- 
ilar to the 43 and 100 GHz lightcurves from the unshocked 
wind emission (see panels h and k in Fig. |9}. 

Radio-to-sub- mm spectra for a variety of viewing orien- 
tations are shown in Fig. 1101 The emission from the WCR 
produces a "bump" on top of the unshocked winds spec- 
trum in a manner which is not too dissimilar to a weak 
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Figure 7. Intensity images from model cwb2 at i = 90°. From top to bottom the images are at a frequency of 15, 43, 100, 250, 
and 1000 GHz. From left to right the phase of the observation increases from 0.0 (conjunction, <f> = 0°), to 0.125 (<f> = 315°), to 0.25 
(quadrature, <f> = 270°), to 0.375 (<f> = 225°). The grey scale is the same as in Fig. [6] 



synchrotron source. Note, however, that this work only con- 
siders the free-free emi ssion from the WCR . Such behaviour 
was previously noted in Pittard et al.l (|2006l ) . It is interesting 
that lContreras et al.l (| 19961 ) purposefully observed a sample 
of O and WR stars at 43 GHz to minimize the effects of 
non-thermal emission on determinations of mass-loss rates. 
In this respect Fig. \TU\ shows that even at 43 GHz the con- 
tribution to the total flux due to thermal emission from the 
WCR may not be negligble, complicating this goal. 

Nevertheless, at high frequencies the thermal emission 
from the unshocked winds always dominates the emission. 
This occurs at v > 100 GHz in model cwb2. Naively, there- 



fore, one might expect that sub-mm observations may be 
preferred to radio observations in respect of determining the 
mass-loss rates of the winds. However, there is then the ad- 
ditional complication that the acceleration regions of the 
winds may be being probed. Another problem is that the 
effects of clumping in the winds may be more severe at mm 
wave l engths than in the r adio (see e.g. lRunacres fc Blommel 
ll99d : lBlomme et al.|[2003 ). Moreover, the flux from sub-mm 
observations is typically not as well calibrated as the flux 
from radio observations. 

The lower panels of Fig. [10] show the spectral indices as 
a function of frequency. a W mds (as shown by the pink line) 
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Figure 8. As Fig. but for an inclination angle i = 30°. From top to bottom the images are at 15, 100, and 1000 GHz, and from left 
to right the orbital phases are 0.0, 0.125, 0.25, and 0.375. The stars and the WCR again rotate anticlockwise, while the grey scale is the 
same as in Fig. \E\ 



is consistently around +1.0, irrespective of the orientation 
of the system to the observer. At lower frequencies a w inds is 
slightly larger, though this may simply reflect some of the 
flux loss which occurs from the model due to the finite grid 
size. At high frequencies the emission from the WCR has a 
spectral index qwcr ~ —0.1. awCR increases at lower fre- 
quencies as the circumstellar environment absorbs some of 
the intrinsic free-free emission from the WCR. The spectral 
index of the total free-free emission, a, lies between a w i n d s 
and awCR- Between 43 and 1000 GHz, a = +0.69, +0.67, 
and +0.65 in Figs. HOf dVffl. These a re comparable to the 
values observed from some stars (e.g. Williams et al. 1990l; 
Leitherer fc Robert|[l99ll ; lAltenhoff et al. Hl994l ; lNugis et all 
19981 ). However, the spectral index from the unshocked 
winds over this frequency range is steeper (a w i n ds = +0.93, 
+0.99, and +1.00), being offset by the flatter spectra from 
the WCR emission («wcr = +0.16, +0.00, and -0.04, in 
Figs. HOf dVff). respectively). 

The variation of the total thermal emission for an ob- 
server in the orbital plane between conjunction and quadra- 
ture is shown in Fig. [5jb). The emission at quadrature is 
dominant over the emission at conjunction over the entire 
frequency range considered. 



4.2.3 Model cwb3 

The left panel of Fig. [11] shows intensity images of the free- 
free emission from model cwb3 (an unequal winds simula- 
tion) at 43 GHz for an observer directly above the orbital 
plane. A comparison against the 43 GHz synthetic image 



in Fig. [6] reveals several differences. Firstly, the position of 
maximum brightness (at the apex of the WCR) is differ- 
ent in the two images: in Fig. 1111 the apex of the WCR is 
pushed closer to the 08V star (which is the southern star 
in this image). Secondly, the emission from the downstream 
arms of the WCR is also in different locations, reflecting 
the changed positions of ram-pressure balance between the 
two winds. Finally, there is a greater emission contrast be- 
tween the leading and trailing edges of the leading arm of the 
WCR, which reflects the changes in the underlying density 
distribution (see Paper I). Conversely, there is also a reduced 
emission contrast between the leading and trailing edges of 
the trailing arm, compared to that from model cwb2 (again 
see Paper I for an explanation of the density distributions 
within the WCRs of these models) . 

The right panel in Fig. [11] shows a synthetic image at 
250 GHz for an observer with i = 30° and = 45°. The 
greater emission from the 06V wind is clearly visible, along 
with the "S"-shaped emission from the WCR. 

Fig. [12] shows intensity images of the free- free emission 
from model cwb3 at 100 GHz as a function of phase for an 
observer located in the orbital plane (i = 90°). Again, sig- 
nificant differences are apparent compared to the 100 GHz 
images from model cwb2 (see Fig. [7]). At conjunction with 
the weaker wind from the 08V star in front (0 = 0°), there 
are two main differences. First, the limb brightened edge on 
the leading arm (to the right side of the image) is projected 
closer to the centre of the image and shows greater curva- 
ture, while the corresponding region of limb brightening on 
the trailing arm is located off the side of the image and is 
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Figure 9. Radio lightcurves of the free-free emission from model cwb2 at 15 GHz (top), 22 GHz, 43 GHz, 100 GHz, 250 GHz, and 1000 GHz 
(bottom) for inclination angles i = 30°, 60°, and 90°. The total free-free emission is shown in the left panels, and the contributions from 
the winds and WCR in the middle and right panels, respectively. In all cases, the observer is located along a direction vector specified 
by 4> = 0°. The stars are at conjunction at phases 0.0 and 0.5, and quadrature at phases 0.25 and 0.75. Note that the y-axes (flux scales) 
arc different in each panel. 



not visible. Second, the reduced radius of the optical depth 
unity surface in the 08V wind creates a smaller silhouette 
against the brighter WCR than occurs in model cwb2. 

Differences also occur at other viewing angles. At phase 
0.125 ((f) = 315°), the WCR is seen mostly face on, with 
the foreground 08V wind silhouetted against it (slightly to 
the right of centre). The vertical curvature of the WCR is 
again clear at phase 0.25 (quadrature), reflecting the wrap- 
ping of the WCR around the 08V star (to the right of cen- 
tre) due to the stronger wind from the 06V star (to the 



left of centre). The apex of the WCR, which occurs closer 
to the 08V star (reflecting the position of the ram pres- 
sure balance of the winds), is seen sideways on. The greater 
emission from the 06V wind, relative to that from the 08V 
wind, is also clearly visible. The double- helix- like structure 
seen from model cwb2 at phase 0.375 disappears in model 
cwb3. Instead, the brightest regions of the projected emis- 
sion from the WCR primarily traces the limb-brightened, 
dense, 08V gas on the trailing edge of the leading arm (see 
Paper I). 
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Figure 10. As Fig. [4] but for model cwb2. The free-free emission from the WCR becomes optically thin at v > 5 GHz (i.e . the f ree-free 
opacity of the shocked plasma has an impact on the spectrum below ~ 5 GHz - see also the discussion in lDouehertv et al.l ||2003| )V 
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Figure 12. Intensity images from model cwb3 at 100 GHz and i = 90°. From top to bottom and left to right the orbital phase is 0.0 
{<j> = 0°, conjunction, 08V star in front), 0.125, 0.25 {<j> = 270°, quadrature), 0.375, 0.5 (<j> = 180°, conjunction, 06V star in front), 
0.625, 0.75 {<j> = 90°, quadrature), and 0.875. The maximum intensity is 1.6 X 10 -10 erg cm - 2 s — 1 Hz - 1 ster - 1 and the major ticks on 
each axis mark out 0.5 mas. 



The images from the second half of the orbit 
(0.5 ^ phased 1.0) resemble those in the upper row of 
Fig. [T2]which show the first half of the orbit. However, some 
small differences are seen, which lead also to differences in 
the lightcurves. For example, the conjunction at phase 0.5 
now has the denser 06V wind in the foreground, whereas 
the weaker 08V wind was in front of the WCR at phase 0.0. 
This changes the nature of the silhouette. 

Fig. [13] shows lightcurves at a range of frequencies from 
radio to sub-mm of the free-free emission from model cwb3. 
The lightcurves are reasonably similar to those shown in 
Fig. [9] from model cwb2, so we show a smaller sample 
of frequencies. The slight differences in the total flux and 



the morphology of the variation have their origin in the 
slighter weaker secondary wind in model cwb3. The max- 
ima now show asymmetrical peaks due to the different wind 
strengths. The highest peak occurs when viewing down the 
trailing arm. This is straighter than the leading arm, as the 
curvature of the WCR opposes the coriolis-induced curva- 
ture, and allows more flux to escape through the low opacity 
plasma within the WCR. 

Fig. [14] shows radio-to-sub-mm spectra from model 
cwb3 for three different viewing orientations. Strong simi- 
larities to the results from model cwb2 are again seen. Since 
model cwb3 consists of 06V and 08V stars, rather than 
two identical 06V stars as in the other models, we also 
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Figure 13. As Fig. [9]but for model cwb3, at v = 15, 100, and 1000 GHz. The theoretical flux from a single, terminal velocity, wind from 
an 06V star is again displayed: the flux from the wind of the 08V companion is 73 per cent of this value. Note that the y-axes (flux 
scales) are different in each panel. 
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Figure 14. As Fig. [4] but for model cwb3. Since the stars differ in this model, the individual flux contributions of the winds from the 
06V and 08V stars are plotted in the top panels. In a) and d) the 08V star is in front. Because of the asymmetry of the WCR the 
spectra at phases 0.25 and 0.75 are different, unlike the case in Fig. [4] 



plot in Fig.[T4]the individual contributions of the unshocked 
winds to the total thermal flux. As expected, the flux from 
the denser 06V wind usually dominates that from the 08V 
wind (see Fig. 114b and c). However, for favourable orienta- 
tions (e.g. when the 08V star is directly in front of the 06V 
star), its contribution to the total thermal flux can be com- 
parable to, and at some frequencies exceed, the contribution 
from the 06V wind (Fig. [Tib'). 



The spectral index as a function of frequency is shown 
in Fig. Q3Jd)-(f). The values of a are very similar to those 
in Fig. UOf dVffl from model cwb2. Note that the high fre- 
quency spectral index from the 08V wind exceeds that from 
the 06V wind. This is because the emission probes deeper 
into the acceleration region of the 08V wind at a given 
(high) frequency. 

Fig. [5jc) shows the variation of the total thermal emis- 
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sion for an observer in the orbital plane between conjunction 
(phase 0.0, 08V star in front) and quadrature (phase 0.75). 
Again, the behaviour is similar to that from model cwb2, 
though the variation between phases is slightly weaker. The 
flux, of course, is also reduced. 

4.2.4 Model cwb4 

Unlike the previous models which all had circular orbits, 
model cwb4 simulates a CWB with an eccentric orbit (e = 
0.36). The intrinsic emission now varies with phase, whereas 
it was constant in the circular orbit models cwbl— cwb3. 
Fig. H5l shows intensity images of the free-free emission from 
model cwb4 at 1000 GHz for an observer directly above the 
orbital plane. The images show striking variations in their 
brightness and morphology as a function of orbital phase, 
reflecting the dramatic changes in the WCR during the orbit 
(see Paper I for full details of the hydrodynamics). 

Discussion of this figure begins, most easily, at orbital 
phase 0.5 (apastron). At this phase we see the emission re- 
gions of the two winds (located very close to the stellar sur- 
faces due to the high frequency) , plus between them emission 
from the hot WCR. There are also small but very bright re- 
gions of emission within and close to the WCR. These are 
very dense, cold, clumps of previously shocked gas which 
formed during the previous periastron passage when the 
WCR became highly radiative. Most of the dense, cold, gas 
which formed during this time has now been cleared out of 
the inner parts of the system. The dense clump and its tail 
visible in the top left of the image are actually located out- 
side of the WCR, since the high inertia of the clumps causes 
them to move on almost ballistic trajectories, while the hot 
gas in the WCR responds more rapidly to the changing po- 
sitions of the stars (see Paper I for more details). 

The process of clearing out the remaining dense clumps 
continues as the stars advance in orbital phase. By phase 
0.7, no clumps are left within the central regions of the sys- 
tem. At this point the stars are moving in a slow waltz to- 
wards each other. The density within the WCR steadily in- 
creases as the stellar separation reduces, resulting in brighter 
emission from the WCR, so that by phase 0.95 the apex of 
the WCR is brighter than the unshocked winds themselves, 
despite the high frequency! The spatial distribution of the 
emission from the WCR also reflects the WCR's increasing 
curvature as periastron is approached. 

By periastron the gas in the densest parts of the WCR 
in model cwb4 has undergone substantial radiative cooling. 
Cool (T ~ 10 4 K), dense gas forms in the WCR between the 
stars and for some distance downstream (see Figs. 11 and 15 
of Paper I). This gas is optically thick to frequencies exceed- 
ing 1000 GHz. The enhancement in the density of this gas 
relative to that of the neighbouring unshocked winds means 
that its emission coefficient is significantly higher, as indi- 
cated by the images shown in Fig. [15] Further downstream, 
the plasma in the WCR remains hotter, with some gas at 
temperatures exceeding 10 7 K (see Fig. 15(e) in Paper I). 
This plasma is optically thin to emission above v > 5 GHz, 
and its emission coefficient is much smaller than that of the 
optically thick cold gas nearer the apex of the WCR. The 
wide variation in plasma density and temperature through- 
out the WCR demonstrates that models of the emission from 



such systems must be based on time-dependent hydrody- 
namical calculations. 

Even more mass within the WCR has cooled back to 
the wind temperature (10 4 K) by phase 0.05. Newly shocked 
gas in the WCR becomes largely adiabatic by phase 0.2, as 
the stars continue to separate and the pre-shock velocity of 
the winds increases. No new clumps form after this phase 
until the next periastron passage. Instead, the slow process 
of clearing the existing clumps out of the system begins. 

It is likely that the majority of the emission between 
orbital phase 0.7 — 1.05 is captured by the numerical grid, 
but it is clear that significant flux losses occur from the 
simulation between phase 0.05—0.3 as optically thick clumps 
leave the finite boundaries of the grid. This effect is likely 
to be greatest at the lower frequencies since the clumps stay 
optically thicker for greater outflow distances, and should be 
kept in mind when examining the remaining figures in this 
section. 

Fig. [TB] displays lightcurves of the free-free emission 
from model cwb4. Since the typical temperature of the 
plasma in the WCR varies between "cold" near periastron, 
and "hot" for the rest of the orbit (excepting the dense 
clumps formed within the WCR), we have chosen to sep- 
arately identify the emission from cold (T < 10 5 K) and hot 
(T ^ 10° K) gas, rather than from the unshocked winds and 
the WCR. At times when the WCR is hot (such as at apas- 
tron), this leads to the same identification of gas as in the 
previous sections. However, during the small range in phases 
around periastron, our analysis groups the cold gas within 
the WCR with the unshocked gas from the winds, rather 
than with the hotter parts of the WCR on the leading arms 
and downstream. Hence, the "cold" and "hot" components 
cannot be directly compared to the "winds" and "WCR" 
components from model cwbl, although as noted above it is 
still reasonable to compare to model cwb2 at other phases. 
Note that Fig. ll6l shows observed fluxes, so, for example, the 
contribution of hot material to the flux includes absorption 
from overlying cold material in the circumstellar environ- 
ment. 

Fig. U6l shows that the emission typically reaches a max- 
imum around periastron, as was already obvious from the 
images at v — 1000 GHz and i = 0°, presented in Fig. 1151 
The varying circumstellar absorption due to the passage of 
the stars in front of the WCR complicates the interpretation 
of the lightcurves when the observer is in the orbital plane. 
Therefore, we shall concentrate on the lightcurve for an ob- 
server at i = 0° . At 43 GHz the emission peaks just prior 
to periastron, while at 1000 GHz the maximum occurs just 
after periastron. There are two reasons for this behaviour. 
One reason is related to the changes with orbital phase of 
the optical depth of sightlines through the WCR, and is dis- 
cussed in more detail below. A minimum centered on phase 
0.1 also occurs in the 43 GHz lightcurve - this is primar- 
ily due to the high optical depth of the cold dense plasma 
formed in the WCR during the periastron passage. Most of 
the emission from this material is therefore absorbed. By 
phase 0.3 many parts of the WCR become optically thin, 
and the emission from background material starts to be- 
come visible. However, the emissivity of this material is also 
dropping, and this becomes the dominant factor affecting 
the observed flux at subsequent phases. The 43 GHz emis- 
sion reaches a minimum around phase 0.5 — 0.6, and then 



Radio emission models of colliding wind binaries 17 



- 1 1 ' 1 1 ■ 

1 0.0 1 

IS: 

- . . i . . . i . . . i . i i iM - 


r _u i | i i i | i i i | i i i | i i t 

■ | 0.05 " 

1 

r ,i . r , , , 1 1 M, 


- i 1 1 1 i - 

1 0.1 

, , 1 , , , 1 , , , 1 , , , 1 , , , 1 , , , 


lo, Wfr* 

- t \ 


, | , , i | ,. , , | , , , | , , , | , 

10.3 fc* 

» , , 1 , , , 1 , , , 1 , M\ , 1 , , - 


pi i | i i i | i i i | i i i | i i i | i i 
-OA Z 

-SM - 


• • ' 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 _ 

_ 0.5 

• ~~~ 

1,1 


,,,,,,,, 

_ 0.6 1 

: • " 

».., .. i « „ . i ,. ,. ... i , „ ... i , .. . t ,. .. z 


_ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 _ 
- 0.7 

; * ; 

-ii.ii- 


; 1 1 i 1 1 1 i 1 1 1 i 1 1 1 i 1 1 1 i 1 1 _ 
10.8 1 


; 1 I 1 1 1 I 1 1 1 I 1 1 1 jf" 1 1 I 1 1 _ 

10.9 1 

\ «P \ 

: ^ 


, 1 , 1 , 1 , 1 , 1 I [ 1 1 , 1 1 , , 1 1 I J 
10.95 1 

; / ; 

. , , i , , , i , , , f- i , , , 



Figure 15. Intensity images from model cwb4 at 1000 GHz for an observer with i = 0°. The orbital phase is noted on each panel. The 
maximum intensity (black in the images) is 8.8 X 10 — 8 erg cm -2 s — 1 Hz -1 ster - 1 . The major ticks on each axis mark out 0.2 mas. 



subsequently increases as the stars approach each other, due 
to the increase in density of the hot plasma in the WCR. 

In contrast, the 1000 GHz lightcurve of the total emis- 
sion (Fig. 1161 ) does not display a minimum at phase 0.1. 
This is partly because some parts of the WCR remain opti- 
cally thin at this phase (although a small reduction in the 
emission from "hot" plasma relative to the trend is never- 
theless visible in Fig. 1161). It is also partly because some 
of the emission at 1000 GHz is from the unshocked stellar 
winds (which is responsible for some of the emission seen in 
Fig. 116k ). and the stars are near quadrature at phase 0.1. 

Fig. [17] shows radio-to-sub-mm spectra from model 
cwb4 as a function of orbital phase for an observer in the 
orbital plane at (f> = 0° . The spectrum from the WCR is op- 
tically thin between phase 0.5 and 0.9 over the displayed fre- 
quency range. Comparing against the spectrum from model 
cwb2 in Fig. [TOTa). we see that both the "hot/ WCR" and 
the "cold/ winds" flux are in good agreement at high fre- 
quencies. However, there is less emission from model cwb4 
than from model cwb2 at low frequencies. At 10 GHz, the 
"winds" flux from model cwb2 is about three times greater 
than the "cold" flux from model cwb4, while the "WCR" 
flux is more than an order of magnitude greater than the 
"hot" flux. This lack of emission from model cwb4 at the 
lower frequencies reflects the smaller grid compared to that 
of model cwb2. 

Since a large fraction of the emission from the WCR 
is optically thick at periastron, it is not surprising to find 
that the spectral index of the "hot" emission at phase 1.0 
becomes positive like the unshocked winds (qwcr ~ +0.75 
between 15 and 100 GHz - see Fig. [T8p). In fact, the "hot" 



spectrum is remarkably straight at periastron, with very lit- 
tle curvature (see Figs. 117b and k, and Fig. 118b ). 

Due to the dramatic cooling at periastron, by phase 0.1 
the majority of the gas within the WCR is now cold. This 
results in a much lower flux from the "hot" component of 
the emission, which is now mostly low density gas on the 
leading edges of the WCR and a narrow strip of rapidly 
cooling gas near the apex and trailing edges of the WCR. 
Conversely, the emission from the "cold" component is sig- 
nificantly bolstered by the substantial amount of gas which 
has cooled within the WCR, being 2.5 x more luminous at 
100 GHz than the emission at periastron. 

The change in viewing angle into the system also causes 
the low frequency emission to drop (see e.g. panel 1 of 
Fig. I17|l . A gradual recovery occurs in the low frequency 
emission up to phase 0.3, after which it declines again. These 
changes are largely driven by changes to the optical depth to 
the "hot" flux, which shows more variability at low frequen- 
cies than the "cold" flux. The level of variability is expected 
to be smaller in calculations where the hydrodynamical grid 
extends to larger radii. 

Although some parts of the WCR remain hot even at 
periastron (for instance the leading edge of the shocks), the 
majority of the mass within the WCR is optically thick at 
phases 0.0 and 0.1. At the highest frequencies considered, the 
WCR is only marginally optically thick at phase 0.1, with 
the optical depth through its surface occasionally exceeding 
unity, but more typically being less than this. The WCR 
becomes increasingly optically thin thereafter, as the stars 
continue their separation. 

Just as the flux from the hot part of the WCR varies 
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Figure 16. Radio lightcurves of the free-free emission from model cwb4 at 43 GHz (top), 100 GHz, 250 GHz, and 1000 GHz (bottom). 
The emission for three different viewing orientations is presented: in the orbital plane (i = 90°) with = 0° (blue) and <f> = 90° (pink), 
and directly above/below the orbital plane (i = 0°, cyan). For the blue curves the stars are at conjunction at phases 0.0 and 0.5, and 
quadrature at phases 0.14 and 0.86. For the pink curves the stars are at conjunction at phases 0.14 and 0.86, and quadrature at phases 0.0 
and 0.5. The total free-free emission is shown in the left panels, and the contributions from the cold and hot gas in the middle and right 
panels, respectively. Note that what is shown in the middle and right panels differs from previous similar figures (see text for details). 
Note also that the y-axes (flux scales) are different in each panel. 



with phase, so does the contribution from the unshocked 
winds plus cold parts of the WCR. At phase 0.5 this part 
of the emission has a straight slope with a co id ~ +1.0 be- 
tween 200 and 1000 GHz. The flux and spectral slope are 
relatively steady up to phase 0.9. At periastron, the first gas 
in the WCR to cool back down to 10 4 K creates a notice- 
able "bump" on the cold gas spectrum at v > 70 GHz (see 
Fig. 118b ) . The extra emission becomes more extended in its 
range of frequencies as more and more gas passing through 
the WCR cools back to the temperature of the unshocked 
winds, and reaches its maximum contribution at phase 0.1, 
with 6x as much flux from cold gas compared to phase 0.5 at 
200 GHz. Thereafter the bump declines as the cold plasma 
in the WCR is ablated and/or expelled out of the system, 
so that the "cold" spectrum is almost entirely generated by 
emission from the unshocked winds by apastron. 

These changes to the spectra also affect the spectral 
index of the emission, as shown in Fig. 1181 with the most 
rapid changes to a occuring around periastron passage. At 
phase 0.9, the spectral index of the total emission increases 
from +0.26 at 55 GHz, to +0.95 at 1 THz, due to the steady 
increase of flux from the cold unshocked winds relative to 
the hot WCR. At periastron, the increasing opacity of the 



hot WCR leads to a significant increase in Qh ot to ~ +0.8 
between 20 and 400 GHz. The enhancement to the emission 
around 100 GHz due to the formation of cold gas increases 
Qcoid to nearly +1.8 at 100 GHz. The spectral index of the 
total emission is in the range +0.9 to +1.3 between 20 and 
600 GHz. By phase 0.1, a monotonically decreases from low 
to high frequencies. However, the computed values of the 
spectral index at low frequencies shown in Fig. I18f c) are 
unlikely to be quantitatively correct, given the known un- 
derestimate of the flux due to the finite grid size. 

Fig. [19] shows the variation of the total thermal flux 
with orbital separation for an observer with i — 0° . Bearing 
in mind that our simulations likely underestimate the true 
flux from this system between phase 0.05 — 0.3, it neverthe- 
less reveals that there is a strong hystersis to the emission at 
all frequencies. It is also interesting that at low frequencies 
the emission is stronger as the stars approach periastron, 
whereas at high frequencies the emission is stronger as the 
stars recede from periastron. This is because at low frequen- 
cies the winds are very opaque and the emission most easily 
escapes when the WCR is full of hot plasma, which is the 
case in the second half of the orbit. In contrast, at high 
frequencies the winds are always transparent, and stronger 
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Figure 18. The spectral index as a function of frequency from the spectra in Fig. ll7l at orbital phase 0.9, 1.0, and 1.1 (which encompasses 
the greatest spectral changes). The observer has a viewing angle i = 90° and <j> = 0°. 
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Figure 19. Hystcrsis of the thermal free-free emission as a function of orbital separation for a variety of frequencies ranging from 15 GHz 
(a) to 1000 GHz (f). The observer is located directly above or below the orbital plane at i = 0°. The phases of some of the points are 
marked. The lower horizontal line shows the theoretical flux calculated from the terminal speed wind of a single 06V star. The upper 
horizontal line shows twice this value. 



intrinsic emission from the cold, dense, plasma which exists 
in the WCR following periastron passage can be observed. 

Fig. [19] also shows the theoretical flux from a single 
06V wind which accelerates instantaneously to its terminal 
velocity (lower horizontal line). The upper horizontal line 
shows twice this value. We have already seen from many of 
the lightcurves in this paper that at the lower frequencies 
the models at times significantly underestimate the true flux 
due to the finite extent of the numerical grid, and we again 
see this here. At higher frequencies (e.g. v > 100 GHz in 
model cwb4) this is no longer an issue, at least over orbital 
phases ^ 0.5 — 1.05. 

The flux at 43 GHz from model cwb4 typically exceeds 
twice the theoretical single star flux, although it is below 
this value at phases 0.05 and 0.1 when the stars orbit deep 
into each other's wind - at these times the optically thick 
contribution from the WCR is not able to offset the loss 
of flux from the unshocked winds due to the presence of 
the WCR (however, particularly at phase 0.1, we expect 
the true flux to be significantly underestimated due to op- 
tically thick, cool, clumps flowing through the grid bound- 
aries). By 250 GHz, the total flux from model cwb4 always 
exceeds twice the expected single star value. This is again 
because the characteristic region of emission is now close 
enough to the stars that the WCR no longer "cuts" out sig- 
nificant parts of this region. At phase 0.5 the extra emission 
from the WCR, which is mostly from hot gas, is negligible 
(see Fig. I17f). However, as the stars approach periastron the 
emission from the WCR becomes significant as the plasma 
within it increases in density and decreases in temperature, 
causing the flux rise shown in Fig. ll9f e). The dip in the flux 
at periastron is caused by the foreground wind occulting the 
apex of the WCR where the brightest emission arises. 

iPittard et ai] (|2006t ) noted that the thermal radio flux 
from the WCR in an adiabatic system scales as cQ p , so that 



we would expect a variation of just over a factor of two from 
our model. However, the variations from model cwb4 are 
obviously much greater. For example, Fig. I16f l) shows that 
for an observer with i — 0° , the flux increases by a factor of 
3.4 between phases 0.6 and 0.9, compared to a d~p increase 
of 1.7. This is because the pre-shock velocities change as the 
stellar separation varies, affec ting the post-shock t empera- 
ture, whereas the analysis in IPittard et al.l (120061 ) is only 
valid in the scale-free limit where the winds collide at con- 
stant velocity. 

Finally, we note that the brightness of the WCR emis- 
sion is highly dependent on its density. In models cwbl 
and cwb4 the cooled postshock gas becomes very dense, 
partly because of the assumption that there is no signifi- 
cant magnetic pressure which would otherwise support this 
gas against collapse. If magnetic pressure were to become 
significant, the densities in the WCR would be limited and 
the emissivity of the WCR would be reduced. 

4-2.5 Clumpy winds 

Figs. [2UJ and [5T] show how the lightcurves and spectra from 
model cwb2 are modified if it is assumed that the unshocked 
winds (and at times also the shocked gas in the WCR) are 
clumpy (note that the mass-loss rates remain the same). 
Fig. 120T b) shows how clumping with a volume filling factor 
/winds = 0.1 increases the flux from the unshocked winds 
compared to the smooth winds case. The flux increases by a 
factor of 2.8, which is a little below the expected increase of 
4.64 from the relation S v oc M 4//3 oc / _2/ ' 3 , and is likely due 
to the greater loss of flux from the numerical grid because 
of its finite size. 

Fig. I20f c) shows how the flux from the WCR varies with 
clumping. Comparing the smooth to the / w i n ds = 0.1 results, 
we see that the larger optical depth unity surfaces in the 
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Figure 20. Comparison of 43 GHz lightcurves from model cwb2 with smooth and clumpy winds for an observer with i = 90° and 

(f> = 0°. The stars are at conjunction at phases 0.0 and 0.5, and quadrature at phases 0.25 and 0.75. Note that while the flux from clumpy 

2/3 1 

unshocked winds increases as / winds , the flux from a clumpy WCR increases as / WCR , since the former is optically thick, while the latter 

is optically thin. Panel a) plots the total flux, panel b) the flux from the unshocked winds, and panel c) the flux from the WCR. Note 

that the y-axes (flux scales) are different in each panel. 
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Figure 21. Free-free radio spectra from model cwb2 as a function of the clumping within the unshocked winds and the WCR, for 
i = 90° and tf> = 0°. a) Comparison of the smooth winds and WCR case (/ w i n( js = /wCR = TO) with clumpy winds but a smooth 
WCR (/ w i n( js = 0.1, /wCR = 1-0). b) As a), but now the WCR is somewhat clumpy too, although the WCR remains much smoother 
than the winds (/ w i n( js = 0.1, /wCR. = 0.5). c) As b), except the clumping within the WCR remains as severe as in the winds 
(/winds = /wCR = 0.1). Note that in the smooth (clumpy) winds case the fluxes are underestimated for v < 20 (50) GHz, because of the 
finite extent of the numerical grid. 



winds of the latter calculation lead to greater attenuation 
of the WCR flux. Hence the observable flux from the WCR 
can be affected by changes to the clumping in the unshocked 
winds, even though in both of these calculations the WCR 
is assumed to be smooth and the intrinsic emission from the 
WCR is identical. 

IPittardl (|2007l ) showed that an adiabatic WCR rapidly 
smoothes out clumps if their density contrast and size 
is not too high, so that the destruction timescale of the 
clump is shorter than the dynamical timescale for material 
within the WCR to flow out of the system. However, this 
work also revealed that the plasma in the WCR contains a 
multitude of weak shocks, and is far from being perfectly 
smooth. To account for this, we have also computed calcu- 
lations where there is some "residual clumping" within the 
WCR. Fig. HJJc) shows that the flux from the WCR grad- 
ually increases as the WCR is made more clumpy. When 
/wcr = 0.5, the observed flux from the WCR at quadra- 
ture is roughly equal to that from the smooth winds case, 
with the extra absorption from the unshocked winds (with 
/winds = 0.1) being offset by extra intrinsic emission from 
the WCR. 

At the extreme limit of there being no "smoothing" ef- 
fect by the WCR (i.e. /wcr = /winds < 1), the observed flux 
from a clumpy WCR can dramatically exceed its smooth 
counterpart, as shown in Fig. I20f c) when /wcr = /winds = 
0.1. However, we regard this possibility as extremely un- 
likely in systems with an adiabatic WCR, and more typically 
would expect /wcr >> /winds when /winds << 1- Indeed, 



the level of the co ntinuum i n the theoretical X-ray spectra 
shown in Fig. 3 of iPittar d (2007) shows that the intrinsic 
flux increases by only 20 per cenlQ equivalent to the WCR 
smoothing the clumping from /winds ~ 0.1 to /wcr ~ 0.8. 

Fig. [21] shows the spectra obtained from calculations 
where the winds are significantly clumped (/winds = 0.1) and 
the degree of clumping within the WCR is varied between 
/wcr = 1-0 (i.e. a smooth WCR) and /wcr = /winds = 0.1 
(i.e. the WCR has no smoothing effect on the clumps). In 
Fig. I21f a) we see that clumping increases the free-free flux 
from the winds, at the expense of increased absorption of 
the free-free flux from the WCR (particularly at the lower 
frequencies). Increasing the clumping within the WCR has 
little noticeable effect at first (see Fig. 121b ). but large en- 
hancements in the flux are apparent when the clumping 
within the WCR is extreme (see Fig. 121b). In fact, in such 
cases the increase in flux from the WCR exceeds the increase 
seen from the winds, due to the differences in optical depth 
within these regions. 

Figs.[2D]and[2Udemonstrate that even with considerable 
clumping in the winds, free-free emission from the WCR can 
still escape and create variability with orbital phase. 



5 ft was this very behaviour which led Pittard (2007) to note that 
"clumping-indcpendent" measurements of the mass-loss could po- 
tentially be obtained from the X-ray flux. 
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5 COMPARISON WITH OBSERVATIONS 
5.1 O-star binaries 

The parameters in model cwbl are similar to those in 
HD 215835 (DHCep), HD 165052, and HD 159176. Cur- 
rently, none of these s ystems is detected i n the radio or 
sub-mm bands, though iBieging et al.l (|l989f ) note 6 cm up- 
per limits of 0.27 mjy for HD 165052, and 1.05 mjy for 
HD 159176. Since we do not capture all of the flux from 
model cwbl (particularly at the lower frequencies investi- 
gated), a direct comparison is not possible. However, we can 
obtain a very conservative upper limit to the 5 GHz flux from 
model cwbl by extrapolating the 1000 GHz flux to 5 GHz 
with a spectral index of +0.6 (the true spectral index may 
be significantly steeper than this). This yields a 5 GHz flux 
of w 0.06 mjy, which is still consistent with the observa- 
tional upper limits noted above, even if moderate clumping 
was added to our models. 

However, these systems, and others like them, should 
be detected by the next generation of interferometers. 
For instance, the EVLA will have a point source sensi- 
tivity better than 1 /xjy between 2 — 40 GHz (la rms in 
ss 9hrs). e-MERLIN's specifications are almost as good. 
The Square Kilometre Array (SKA), which will operate 
over 0.1 — 25 GHz, will probe even deeper. ALMA's band 3 
(84— 116 GHz) also has a continuum sensitivity of 0.06 mjy, 
and should be just sensitive enough to detect these systems. 
Obviously the chances for detection improve if the winds are 
clumped. 

Model cwb2 is similar to HD 93161A, and also 
shares some similariti e s with Plaskett's star (HD 47129). 
iBenaglia fc Koribalskil (120041 1 failed to detect HD 93161 A 
with ATCA at 3 and 6 cm. Observations centered on 
HD93129A (about 3 arcmin from HD93161A) reveal that at 
the position of HD 93161A the rms noise is ^ 0.2 mjy/beam 
at 3 cm, and ^ 0.25 mjy/beam at 6 cm, giving 3cr upper 
limits of 0.6 mjy and 0.75 mjy, respectively. These limits 
are consistent with the fluxes predicted from model cwb2. 
On the other han d, Plaskett's star was detected at 5 GHz by 
IPersi et al.l (119881 ). who reported a flux of 0.2 mjy. This flux 
is higher than we see from model cwb2, and may indicate 
significant clumping within the winds of this system. The 
fact that Plaskett's star is at a distance of 1.5 kpc (rather 
than the 1 kpc of the models in this work) makes this differ- 
ence even worse. 

Model cwb4 bears some similarities to some well-known 
O+O systems with eccentric orbits, including HD 152248, 
HD 93205, HD 93403, CygOB2#8A, and tOri. Of these sys- 
tems, HD 152248, HD 93205, and HD 93403 are currently not 
detected in the radio. HD 152248 fell within the field of view 
of A TCA observations of Sco OBI (|Setia Gunawan et al.l 
2003), and 3cr upper limits (at both 3 cm and 6 cm) are 
0.24 mjy (Stevens, priv. communication, 2009). In contrast, 
Cyg OB2 #8A is a bright, non-therm al, and strongly vari- 
able radio source |Bieging et al.ll 19891). Following the discov- 
ery o f this system's b inarity (|De Becker. Rauw fc Manfroidl 
l2004h ,lBl ommel (|2005h showed that the radio flux density un- 
dergoes phase-locked variation^]. The non-thermal emission 
from our models will be investigated in a future paper. 



Other O+O systems with strongly varying non-thermal ra- 



Finally, t Orionis (HD 37043) was detected by Howarth 
& B rown (1991) at 3.6 cm, wit h a flux of 0.046 ± 0.015 mjy 
(see lLamers fc Leithererlll993l ). Since it was observed only 
at one frequency it is not known if this emission is pre- 
dominantly thermal or non-thermal. Moving the system to 
a distance of lkpc yields a flux of ss 0.01 mjy, comparable 
to the predicted 8 GHz flux from model cwb2, and similar 
to the flux from model cwb4 at certain phases. 



5.2 WR and LBV systems 

Although this work has focussed on O+O systems, 
thermal radio emission is also s een from WR and 
LBV stars (e.g. lAbbott et all 1 19861 ; iLeitherer et al.1 Il997l ; 
iDoughertv fc Clark|[20o3 ). However, there are few accounts 
in the literature on the variability of such sourcefl For in- 
stance, despite a large multi-frequency monitoring campaign 
of 7Velorum, little evidence for variation was found (Sean 
Dougherty, priv. communication, 2009). 

Radio variability at the 10 — 30 per cent level, has, 
however, been recently reported from four WN9 stars in 
the Arches clus ter (sources AR1, AR3, AR4, and AR8: 
lLang et al.ll2005T ). All but AR1 have a positive spectral in- 
dex (a > +0.6; observational uncertainities do not allow 
an accurat e meas urement of the spectral index of AR1). 
lLang et all (j2005r ) suggested the variability might indicate 
the presence of a time-variable non-thermal component, 
and/or changes to the underlying mass-loss rate or veloc- 
ity of the wind. However, as this work demonstrates, an al- 
ternative explanation is that these sources are binaries, and 
that the variability is caused by changes to the circumstellar 
opacity along sightlines into the system as the stars orbit 
each other, or due to a non-zero orbital eccentricity. This 
view is supported by the fa ct that AR1 and AR4 (iden tified 
as sources A1N and A1S bv lWang. Dong fc Lanell2006l ) have 



strong and hard X-ray emission with a prominent FeK line 
at 6.7 keV, all of which are characteristics of colliding wind 
binarie^f). 

Three point-like radio sources were a lso detected in th e 
Quintuplet cluster (QR6, QR8, and QR9: lLang et al.ll2005l V 
While these all have positive spectral indices indicative of 
a strong stellar wind, we have demonstrated in this work 
that the bulk of the emission could actually arise from a 
radiative, and optically thick, WCR. In Westerlundl, six 
WR stars are detected in the radio, but only one has a ther- 
mal spectrum (WRL; Simon Clark, priv. communication, 
2009). Note, however, that WR L looks like it mi ght be a 
CWB from it's X-ray properties (|Clark et al.ll200Sl ) . 

Variable thermal emission with a > +0.6 has also been 
detec ted from WR89, WR113, and WR138 (|Montes et al.1 
2009). For WR89, the spectral index of the emission, 
a > +0.6, while for WR113, a > +0.9, and a > 



dio emission are HP 168112 llBlomme et al.l |2005|), HP 167 971 
llBlomme et al.ll2007t) , and CygOB2#9 l lVan Loo et al.ll2008l) . 

7 Of course, the variability of many non-thermal sources is much 
better studied. 

8 On the other hand, source AR6 clearly shows a non-therma l 
radio spectrum (lLang. Goss fc Rodri'guej200ll;lLang" et al .120051), 
and s ince it is associated with strong X-ray emission dWang et al.l 
120061) . it is also likely to be a colliding wind binary. 



Radio emission models of colliding wind binaries 23 



+0.7 for WR138. WR8 9 is classified as a visual bi- 
nary (|van der Huchj|200ll ). while WR113 i s a WC8+Q8- 
9 sys tem with an orbital period of 29.7 d ()Niemela et alJ 
1 19991 ). WR 138 is a s pectroscopic binary with a period of 
1538 d jAnnuldll99(ih . though this system may also har- 
bour a close companion with a short period o f order days 
|Lamontagne et al.lll982l ; iMoffat fc Sharalll98rt) . It is likely 
that the obser ved variability is rel ated to the binarity of 
these systems. iMontes et al.l (| 2009? ) also note that the ra- 
dio emission from WR 98a and WR 104 (two of the so-called 
"pinwheel" systems) is variable and the spectral index is 
positive, but less than the value expected from a single wind 
(i.e. a < +0.6). These systems could harbour non-thermal 
emission (i.e. be "composite" systems), or the emission could 
be purely thermal, combining optically thick emission from 
the winds with optically thin emission from the WCR. 

It is also interesting to note that the broad maxima seen 
in the lightcurves of model cwbl are simil ar to those seen in 
the r adio lightcurve of the LBV r/ Carinae (|Duncan fc White! 
2003). The latter is usually explained in terms of time- 
variable ionization of dense gas w ithin the spoked equa torial 
disc seen in the optical (see also lKashi fc So kcr 2007). The 
spatial extent of the brightest part of the radio emission 
spans a width of 1 to several arcseconds, depending on the 
orbital phase. At a distance of 2.3 kpc, this corresponds to 
thousands of AU. However, 3-dimensional simulations of the 
structure of the WCR reveals density enhancements o n simi- 
lar spatial scales (|Parkin et al.ll2009l ; [Gull et al Jl2009h ■ Thus, 
it is possible that at least some of the variability could be 
due to time variations in the properties and orientation of 
the optically thick shocked wind of the primary LBV star. 
Furthermore, ?7Car showed a sharp peak in the 7 mm flux 
density in 2003.5 (and again in 2009.0, Abraham, private 
communication, 2009) . coinciding with periastron passage 
(| Abraham et alj|2005ah . This spike has been interpreted as 
the shocked plasma from the secon dary wind becoming op - 
tically thick at high temperatures (I Abraham et al.ll2005bh . 
It is interesting that we also see a sharp maximum in the 
flux from model cwb4 near periastron. However, for this we 
have a different interpretation. At low frequencies the spike 
in model cwb4 occurs before periastron, when the density 
and thus emissivity of the WCR is increasing. However, the 
WCR is optically thin at this time. In contrast, at high fre- 
quencies (y > 250 GHz), the greatest flux from model cwb4 
occurs after periastron. At such frequencies, the majority of 
the unshocked winds is optically thin, allowing the bright 
emission from the cold dense clumps created during perias- 
tron passage to be seen. It is clear that further modelling 
of n Car is necessary, in order to better understand the ob- 
served flux variations, on both short and long timescales, 
from this fascinating system. 



6 SUMMARY AND CONCLUSIONS 

This work examines the thermal radio-to-sub-mm emission 
from colliding wind binaries, focussing in this first instance 
on short-period O+O binaries (the potential non-thermal 
emission will be studied in due course). We investigate how 
the thermal emission is influenced by the presence of a 
nearby companion star, and the extent to which variabil- 
ity with orbital phase occurs. The emission is calculated 



from grids of density and temperature ge nerated f rom th e 
3D hydrodynamical models presented in IPittardl (1200981 ). 
The models span a variety of different physical conditions, 
including radiative and adiabatic wind-wind collision regions 
in systems with circular orbits, and also a system with an 
eccentric orbit where the WCR repeatedly cycles from ra- 
diative to adiabatic and back to radiative again. 

We find distinct differences in the spectra and 
lightcurves from systems with radiative and adiabatic 
WCR's. In the former (model cwbl), the WCR is optically 
thick to wavelengths as short as sub-mm, and radiates as a 
dense sheet of emission. The WCR dominates the total emis- 
sion from the system, and only begins to show signs of be- 
coming optically thin at v > 1000 GHz. The flux (which has 
a positive spectral index) can be over an order of magnitude 
greater than from the wind of a single star, and highlights 
potentially serious consequences for observational estimates 
of mass-loss rates where sources are assumed to be single 
stars, but are in fact binaries. The flux is highest just before 
quadrature, when the inner regions of the WCR are viewed 
almost face on (due to the large aberration of the WCR), 
and have the largest projected extent. The lightcurve max- 
ima are broad, with minima occuring slightly before con- 
junction. For an observer in the orbital plane, the total flux 
(from the unshocked winds plus the WCR) varies by a fac- 
tor of two or so. Although there is some loss of flux from 
the models, the spectral index of the emission is likely to 
be steeper than +0.6, and varies with orbital phase and the 
orientation of the observer. Synthetic images for an observer 
directly above the orbital plane display an "S" -shaped re- 
gion of emission which traces the location of the WCR. 

In contrast, the emission from adiabatic WCR's is op- 
tically thin, and most easily escapes the system when sight- 
lines to the observer are themselves through the WCR. 
In model cwb2 the WCR dominates the emission at v < 
50 — 100 GHz, with free-free absorption creating a low fre- 
quency turnover typically at v ~ 10 GHz. The emission 
from the WCR typically contributes only 10 per cent of 
the total flux at 1000 GHz (0.3 mm). The lightcurves show 
relatively narrow maxima at phases just after quadrature. 
The broad, and relatively flat minima inbetween the max- 
ima indicate phases where there is considerable absorption 
of the emission from the WCR. There is a gradual drop 
in the degree of variability with frequency as the circum- 
stellar environment becomes less opaque and increasingly 
transparent. At 1000 GHz the vast majority of the WCR 
emission escapes the system. Between 100 and 250 GHz, the 
lightcurves change morphology. At v > 250 GHz, the max- 
ima become almost completely flat-topped, and the minima 
are now caused by occultation of a narrow emission region 
near the surface of each star. The spectral index is steep- 
est at low frequencies (where absorption of emission from 
the WCR is important) and at high frequencies (where the 
emission is almost entirely from deep within the accelera- 
tion region of the individual winds). Images of the emission 
for an observer in the orbital plane reveal an intertwined 
"double-helix" showing limb-brightened parts of the WCR, 
against which the foreground wind is sometimes silhouetted. 
At other viewing angles (e.g. i = 30°) an "S"-shaped region 
of emission is again seen. 

Adiabatic systems where one of the winds is weaker 
than the other generally display similar properties, the 
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main differences being slight asymmetries to the lightcurves 
(e.g. differing heights to the maxima), and a change 
from a "double-helix" emission structure to a single limb- 
brightened structure emitted by the relatively dense trailing 
edge of the trailing arm of the WCR. 

One of the more surprising findings in this work is the 
strength of the flux variations with orbital phase from a sys- 
tem with an eccentric orbit (model cwb4), which exceed an 
order of magnitude from the WCR), and may be almost as 
strong from the system as a whole. This is caused by dra- 
matic changes to the physical properties of the WCR, which 
is highly radiative at periastron, but adiabatic at apastron, 
and changes from optically thick to optically thin, respec- 
tively. This creates a strong frequency-dependent hystersis 
to the emission around the orbit. At low frequencies the 
emission is stronger as the stars approach periastron, due 
to the low opacity through the hot WCR at this time. The 
flux variations easil y exceed the eQp scaling proposed by 
IPittard et al.l (|2006l ), which is invalidated by changes to the 
pre-shock wind speeds as the stars orbit deep within the 
acceleration zone of each other's wind. In contrast, at high 
frequencies the emission is stronger as the stars recede from 
periastron, because of the high intrinsic emission from dense, 
cold, post-shock gas created during periastron passage. Even 
larger variations are seen when the observer is in the orbital 
plane, indicating that changes in the circumstellar absorp- 
tion also play an important role. 

An examination of the effects of clumping reveal that in 
addition to higher fluxes, clumping can also effect the vari- 
ability of radio lightcurves. In the case of strong clumping 
in the unshocked winds but a smooth WCR, the variabil- 
ity decreases. However, if the clumps are not smoothed by 
the WCR (a case which we consider rather unlikely), the 
variations of the flux with orbital phase can be enhanced 
compared to the smooth winds case. 

We find that the fluxes from our models are consis- 
tent with the mostly upper limits which exist to date from 
the O+O systems which they most resemble. On the other 
hand, there are many instances where radio variability is ob- 
served from WR sources with positive spectral indices, some 
with a > +0.6. Many of these sources are likely to be bina- 
ries, with the variation linked to changes in the absorption 
through the circumstellar material, as the orientation of the 
system to the observer changes, or perhaps due to changes 
in the separation of the stars. However, we also caution that 
systems with positive spectral indices are not necessarily 
thermal - instead, one could be witnessing the absorption 
of a non-thermal component. This possibility must be care- 
fully considered, especially if the spectral index is calculated 
at relatively low frequencies (e.g. between 5 and 8 GHz). To 
safely distinguish between thermal and non-thermal emis- 
sion in close, spatially unresolved systems where the spectral 
index is positive may require flux measurements over a wide 
frequency range, additional observations at X-ray and 7-ray 
energies, and detailed comparisons with theoretical models. 

It is difficult to predict how the fluxes from the mod- 
els investigated in this work will scale to other systems (e.g. 
WR+O systems) , since there are so many key variables, such 
as whether the WCR is optically thick or thin, the recent 
thermal history of the plasma in the WCR, the relative sizes 
of the radio photospheres in each wind to the stellar separa- 
tion, the orientation of the system, etc. In many ways this 



work, in revealing some of these complications, demonstrates 
that it may not always be valid to apply so me of the scal- 
ing law s noted in iDouehertv et "all (|2003h and lPittard et all 
(2006). In this respect there is a clear need for further nu- 
merical simulations of CWBs with a range of key parameters 
(mass-loss rates, stellar separations, etc.). Such simulations 
will aid the study and interpretation of the resulting emis- 
sion from these fascinating systems. 

The future commissioning of more sensitive telescopes, 
such as the EVLA, e-MERLIN, ALMA, and SKA, should 
result in detections of short-period O+O systems, thus al- 
lowing the predictions of this work to be tested. Although 
the synthetic images shown in this work are of an angular 
scale which is too small to be resolved with EVLA or ALMA, 
planned upgrades to the VLBA offer the exciting prospect 
of spatially resolving the emission from tight CWBs. The 
EVLA, ALMA, and other facilities may also resolve some of 
the wider O+O (and WR+O) CWBs. 
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